In [ ]:
import matplotlib.pyplot as plt
from matplotlib.colors import LinearSegmentedColormap
import numpy as np
import pyvista as pv
# Définir les couleurs du bas au haut de la colormap de l'image
colors = [
(75/255, 0, 130/255), # indigo
(0, 0, 255/255), # bleu
(0, 255/255, 0), # vert
(255/255, 255/255, 0), # jaune
(255/255, 127/255, 0), # orange
(255/255, 0 , 0), # rouge
(148/255, 0, 211/255), # violet
(1, 1, 1) # blanc
]
cmap_name = 'custom_gradient'
# Créer la colormap
custom_cmap = LinearSegmentedColormap.from_list(cmap_name, colors, N=256)
# Afficher la colormap dans Matplotlib
plt.imshow(np.linspace(0, 1, 256)[None, :], aspect='auto', cmap=custom_cmap)
plt.colorbar()
plt.show()
In [ ]:
import os
import feelpp
from feelpp.toolboxes.cfpdes import *
from tools.lap2D_pinns import Run_laplacian2D, Poisson_2D, PoissonDisk2D
from scimba.equations import domain
class Poisson:
"""
Solves the problem
-Laplacian u = f in Omega
u = g in boundary
- with f,g are set by the user
"""
def __init__(self, dim=2, order=1):
self.dim = dim
self.model = dict()
self.order = order
def genCube(self, filename, h=0.1):
"""
Generate a cube geometry following the dimension self.dim
"""
geo="""SetFactory("OpenCASCADE");
h={};
dim={};
""".format(h, self.dim)
if self.dim==2 :
geo+="""
Rectangle(1) = {0, 0, 0, 1, 1, 0};
Characteristic Length{ PointsOf{ Surface{1}; } } = h;
Physical Curve("Gamma_D") = {1,2,3,4};
Physical Surface("Omega") = {1};
"""
elif self.dim==3 :
geo+="""
Box(1) = {0, 0, 0, 1, 1, 1};
Characteristic Length{ PointsOf{ Volume{1}; } } = h;
Physical Surface("Gamma_D") = {1,2,3,4,5,6};
Physical Volume("Omega") = {1};
"""
with open(filename, 'w') as f:
f.write(geo)
def __call__(self,
h, # mesh size
order=1, # polynomial order
name='Potential', # name of the variable u
rhs='8*pi*pi*sin(2*pi*x)*sin(2*pi*y)', # right hand side
diff='{1,0,0,1}', # diffusion matrix
g='0',
geofile=None,
plot=None,
solver='feelpp'):
"""
Solves the problem where :
- h is the mesh size
- order the polynomial order
- rhs is the expression of the right-hand side f(x,y)
"""
self.pb = cfpdes(dim=self.dim, keyword=f"cfpdes-{self.dim}d-p{self.order}")
self.model = {
"Name": "Laplacian",
"ShortName": "Laplacian",
"Models":
{
f"cfpdes-{self.dim}d-p{self.order}":
{
"equations":"poisson_eq"
},
"poisson_eq":{
"setup":{
"unknown":{
"basis":f"Pch{order}",
"name":f"{name}",
"symbol":"u"
},
"coefficients":{
"c": f"{diff}:x:y" if self.dim == 2 else f"{diff}:x:y:z",
"f": f"{rhs}:x:y" if self.dim == 2 else f"{rhs}:x:y:z"
}
}
}
},
"Materials":
{
"Omega":
{
"markers":["Omega"]
}
},
"BoundaryConditions":
{
"poisson_eq":
{
"Dirichlet":
{
"g":
{
"markers":["Gamma_D"],
"expr":f"{g}:x:y"
}
}
}
},
"PostProcess":
{
f"cfpdes-{self.dim}d-p{self.order}":
{
"Exports":
{
"fields":["all"],
"expr":{
"rhs" : f"{rhs}:x:y" if self.dim == 2 else f"{rhs}:x:y:z"
}
}
}
}
}
fn = None
if geofile is None:
fn = f'omega-{self.dim}.geo'
self.genCube(fn, h)
else:
fn = geofile
##________________________
if solver == 'feelpp':
print(f"Solving the laplacian problem for hsize = {h}...")
feelpp_mesh = feelpp.load(feelpp.mesh(dim=self.dim, realdim=self.dim), fn, h)
self.pb.setMesh(feelpp_mesh)
self.pb.setModelProperties(self.model)
self.pb.init(buildModelAlgebraicFactory=True)
self.pb.printAndSaveInfo()
self.pb.solve()
self.pb.exportResults()
#measures = self.pb.postProcessMeasures().values()
"""
try:
import pandas as pd
df=pd.DataFrame([measures])
print(df)
except ImportError:
print("cannot import pandas, no problem it was just a test")
"""
##______________________
#mesh = self.pb.mesh()
#function_space = self.pb.functionSpace()
#dofs = function_space.dofs()
##________________________
#return mesh, dofs
##________________________
elif solver == 'scimba':
print("Solving using Scimba")
# Define a disk domain
if geofile == 'geo/disk.geo' :
xdomain = domain.SpaceDomain(2, domain.DiskBasedDomain(2, center=[0.0, 0.0], radius=1.0))
pde_disk = PoissonDisk2D(xdomain, rhs= rhs, g= g)
Run_laplacian2D(pde_disk)
# Define a square domain
elif geofile == None:
xdomain = domain.SpaceDomain(2, domain.SquareDomain(2, [[0.0, 1.0], [0.0, 1.0]]))
pde = Poisson_2D(xdomain, rhs= rhs, g= g)
Run_laplacian2D(pde)
##________________________
# Plots
if plot != None:
from xvfbwrapper import Xvfb
import pyvista as pv
import matplotlib.pyplot as plt
#vdisplay = Xvfb()
#vdisplay.start()
pv.set_jupyter_backend('static')
pv.start_xvfb()
mesh = pv.get_reader(f"cfpdes-{self.dim}d-p{self.order}.exports/Export.case").read()
pl = pv.Plotter(shape=(1,2))
pl.add_title(f'Solution P{order}', font_size=18)
pl.add_mesh(mesh[0], scalars = f"cfpdes.poisson_eq.{name}", cmap=custom_cmap)
pl.subplot(0,1)
pl.add_title('f=' + rhs, font_size=18)
pl.add_mesh(mesh[0].copy(), scalars = 'cfpdes.expr.rhs', cmap=custom_cmap)
pl.link_views()
if self.dim ==3:
pl.view_isometric()
else:
pl.view_xy() # if mesh_2D is on the xy plane.
pl.show()
pl.screenshot(plot)
Using device: cpu torch loaded; device is cpu
In [ ]:
import sys
import feelpp
import feelpp.toolboxes.core as tb
# mandatory things
sys.argv = ["feelpp_app"]
e = feelpp.Environment(sys.argv,
opts=tb.toolboxes_options("coefficient-form-pdes", "cfpdes"),
config=feelpp.localRepository('feelpp_cfpde'))
# ------------------------------------------------------------------------- #
# Poisson problem
# - div (diff * grad (u)) = f in Omega
# u = g in Gamma_D
# Omega = domain, either cube or ball
# Approx = lagrange Pk of order order
# mesh of size h
[ Starting Feel++ ] application feelpp_app version 0.1 date 2024-May-21 . feelpp_app files are stored in /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 .. logfiles :/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/logs
In [ ]:
# 2D with varying RHS
P = Poisson(dim = 2)
P(h=0.08, rhs='-1.0-1*y*x+y*y', g='0', order=1, geofile='geo/disk.geo', plot='2d.png')
P(h=0.08, rhs='-1.0-1*y*x+y*y', g='0', order=1, geofile='geo/disk.geo', solver='scimba')
P(h=0.1, rhs='-1.0-2*y*x+y*y', g='0', order=1, plot='f2.png')
P(h=0.1, rhs='-1.0-2*y*x+y*y', g='0', order=1, solver ='scimba')
Solving the laplacian problem for hsize = 0.08...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/geo/disk.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/geo/disk.geo 0.08
+--------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+--------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-17 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 1.012276e-01 | | | |
| | | | h_max : 1.349240e-01 | | | |
| | | | h_min : 8.067356e-02 | | | |
| | | | n_elements : 757 | | | |
| | | | n_faces : 1167 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 411 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +----------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +----------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +--------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +--------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+--------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+==============+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1,0,0,1} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+--------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | -1.0-x*y+y^2 | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+--------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +--------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-17 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 1.012276e-01 | | | | |
| | | | | h_max : 1.349240e-01 | | | | |
| | | | | h_min : 8.067356e-02 | | | | |
| | | | | n_elements : 757 | | | | |
| | | | | n_faces : 1167 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 411 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : 0 | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-17 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 411 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +-------------------------------------------------+ | |
| | | Fields | | |
| | +-------------------------------------------------+ | |
| | | +---------------------------------------------+ | | |
| | | | Potential | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-16 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-16 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-16 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | +-------------------------------------------------+ | |
| +----------------------------------------------------------------------------------------------------------------------------------------------+ |
+--------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 6.73e-01 epoch 0: best loss = 6.73e-01 epoch 1: best loss = 6.34e-01 epoch 2: best loss = 5.72e-01 epoch 3: best loss = 4.43e-01 epoch 4: best loss = 2.83e-01 epoch 5: best loss = 1.39e-01 epoch 12: best loss = 1.33e-01 epoch 13: best loss = 1.33e-01 epoch 14: best loss = 1.25e-01 epoch 15: best loss = 1.13e-01 epoch 16: best loss = 1.05e-01 epoch 17: best loss = 9.61e-02 epoch 18: best loss = 9.10e-02 epoch 19: best loss = 8.59e-02 epoch 20: best loss = 8.46e-02 epoch 22: best loss = 7.89e-02 epoch 23: best loss = 6.18e-02 epoch 24: best loss = 4.50e-02 epoch 25: best loss = 3.54e-02 epoch 26: best loss = 3.49e-02 epoch 27: best loss = 2.60e-02 epoch 28: best loss = 1.25e-02 epoch 34: best loss = 1.23e-02 epoch 36: best loss = 1.06e-02 epoch 37: best loss = 8.09e-03 epoch 44: best loss = 7.14e-03 epoch 46: best loss = 7.08e-03 epoch 47: best loss = 5.63e-03 epoch 48: best loss = 3.78e-03 epoch 55: best loss = 3.77e-03 epoch 57: best loss = 3.29e-03 epoch 58: best loss = 2.72e-03 epoch 64: best loss = 2.68e-03 epoch 65: best loss = 2.26e-03 epoch 68: best loss = 2.19e-03 epoch 71: best loss = 1.92e-03 epoch 72: best loss = 1.86e-03 epoch 74: best loss = 1.67e-03 epoch 77: best loss = 1.58e-03 epoch 78: best loss = 1.54e-03 epoch 80: best loss = 1.39e-03 epoch 81: best loss = 1.30e-03 epoch 83: best loss = 1.26e-03 epoch 84: best loss = 1.25e-03 epoch 85: best loss = 1.24e-03 epoch 86: best loss = 1.15e-03 epoch 87: best loss = 1.05e-03 epoch 88: best loss = 9.94e-04 epoch 89: best loss = 9.66e-04 epoch 90: best loss = 9.14e-04 epoch 91: best loss = 8.89e-04 epoch 93: best loss = 8.03e-04 epoch 94: best loss = 7.97e-04 epoch 95: best loss = 7.17e-04 epoch 97: best loss = 6.51e-04 epoch 98: best loss = 6.27e-04 epoch 99: best loss = 6.14e-04 epoch 100: best loss = 6.02e-04 epoch 101: best loss = 5.20e-04 epoch 102: best loss = 5.06e-04 epoch 104: best loss = 4.72e-04 epoch 105: best loss = 4.23e-04 epoch 106: best loss = 4.13e-04 epoch 108: best loss = 3.83e-04 epoch 111: best loss = 3.53e-04 epoch 112: best loss = 3.50e-04 epoch 114: best loss = 3.12e-04 epoch 116: best loss = 3.10e-04 epoch 117: best loss = 3.00e-04 epoch 121: best loss = 2.80e-04 epoch 124: best loss = 2.60e-04 epoch 126: best loss = 2.53e-04 epoch 127: best loss = 2.45e-04 epoch 130: best loss = 2.45e-04 epoch 131: best loss = 2.38e-04 epoch 132: best loss = 2.35e-04 epoch 135: best loss = 2.16e-04 epoch 140: best loss = 2.15e-04 epoch 142: best loss = 2.08e-04 epoch 143: best loss = 1.98e-04 epoch 147: best loss = 1.97e-04 epoch 149: best loss = 1.91e-04 epoch 151: best loss = 1.81e-04 epoch 152: best loss = 1.75e-04 epoch 156: best loss = 1.68e-04 epoch 166: best loss = 1.67e-04 epoch 167: best loss = 1.54e-04 epoch 177: best loss = 1.50e-04 epoch 182: best loss = 1.47e-04 epoch 183: best loss = 1.45e-04 epoch 185: best loss = 1.45e-04 epoch 186: best loss = 1.44e-04 epoch 188: best loss = 1.41e-04 epoch 189: best loss = 1.39e-04 epoch 212: best loss = 1.39e-04 epoch 220: best loss = 1.32e-04 epoch 227: best loss = 1.30e-04 epoch 228: best loss = 1.25e-04 epoch 236: best loss = 1.15e-04 epoch 284: best loss = 1.14e-04 epoch 292: best loss = 1.12e-04 epoch 293: best loss = 1.05e-04 epoch 297: best loss = 1.02e-04 epoch 302: best loss = 1.00e-04 epoch 312: best loss = 9.58e-05 epoch 317: best loss = 9.49e-05 epoch 323: best loss = 9.45e-05 epoch 329: best loss = 9.11e-05 epoch 335: best loss = 9.02e-05 epoch 349: best loss = 8.95e-05 epoch 358: best loss = 8.39e-05 epoch 368: best loss = 7.76e-05 epoch 496: best loss = 7.62e-05 epoch 500: current loss = 8.00e-05 epoch 511: best loss = 7.46e-05 epoch 512: best loss = 7.42e-05 epoch 527: best loss = 7.27e-05 epoch 541: best loss = 7.24e-05 epoch 552: best loss = 6.89e-05 epoch 573: best loss = 6.81e-05 epoch 574: best loss = 6.77e-05 epoch 593: best loss = 6.58e-05 epoch 595: best loss = 6.55e-05 epoch 599: current loss = 6.67e-05 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
Solving the laplacian problem for hsize = 0.1...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo 0.1
+----------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+----------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-19 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 1.018751e-01 | | | |
| | | | h_max : 1.168628e-01 | | | |
| | | | h_min : 8.366122e-02 | | | |
| | | | n_elements : 248 | | | |
| | | | n_faces : 392 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 145 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+================+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1,0,0,1} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | -1.0+y^2-2*x*y | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-19 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 1.018751e-01 | | | | |
| | | | | h_max : 1.168628e-01 | | | | |
| | | | | h_min : 8.366122e-02 | | | | |
| | | | | n_elements : 248 | | | | |
| | | | | n_faces : 392 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 145 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : 0 | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-19 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 145 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +-------------------------------------------------+ | |
| | | Fields | | |
| | +-------------------------------------------------+ | |
| | | +---------------------------------------------+ | | |
| | | | Potential | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-18 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-18 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-18 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | +-------------------------------------------------+ | |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
+----------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 1.22e+00 epoch 0: best loss = 1.22e+00 epoch 1: best loss = 8.67e-01 epoch 2: best loss = 4.85e-01 epoch 3: best loss = 2.30e-01 epoch 4: best loss = 2.23e-01 epoch 7: best loss = 2.09e-01 epoch 8: best loss = 2.00e-01 epoch 11: best loss = 1.89e-01 epoch 21: best loss = 1.88e-01 epoch 22: best loss = 1.76e-01 epoch 25: best loss = 1.68e-01 epoch 26: best loss = 1.57e-01 epoch 27: best loss = 1.52e-01 epoch 29: best loss = 1.38e-01 epoch 31: best loss = 1.28e-01 epoch 32: best loss = 1.22e-01 epoch 33: best loss = 1.11e-01 epoch 36: best loss = 1.03e-01 epoch 38: best loss = 9.47e-02 epoch 40: best loss = 8.34e-02 epoch 43: best loss = 7.58e-02 epoch 45: best loss = 6.99e-02 epoch 46: best loss = 6.99e-02 epoch 47: best loss = 6.79e-02 epoch 49: best loss = 6.70e-02 epoch 50: best loss = 5.78e-02 epoch 52: best loss = 5.25e-02 epoch 53: best loss = 5.22e-02 epoch 54: best loss = 4.22e-02 epoch 56: best loss = 4.01e-02 epoch 59: best loss = 3.75e-02 epoch 60: best loss = 3.21e-02 epoch 62: best loss = 3.10e-02 epoch 63: best loss = 2.96e-02 epoch 64: best loss = 2.74e-02 epoch 65: best loss = 2.71e-02 epoch 66: best loss = 2.39e-02 epoch 71: best loss = 2.26e-02 epoch 72: best loss = 2.02e-02 epoch 73: best loss = 2.01e-02 epoch 74: best loss = 1.78e-02 epoch 79: best loss = 1.68e-02 epoch 83: best loss = 1.47e-02 epoch 86: best loss = 1.29e-02 epoch 112: best loss = 1.12e-02 epoch 114: best loss = 1.08e-02 epoch 120: best loss = 8.30e-03 epoch 121: best loss = 7.94e-03 epoch 136: best loss = 6.45e-03 epoch 153: best loss = 6.35e-03 epoch 161: best loss = 6.10e-03 epoch 163: best loss = 5.31e-03 epoch 188: best loss = 4.75e-03 epoch 189: best loss = 4.07e-03 epoch 226: best loss = 3.38e-03 epoch 254: best loss = 3.22e-03 epoch 261: best loss = 2.69e-03 epoch 500: current loss = 7.77e-03 epoch 563: best loss = 2.33e-03 epoch 599: current loss = 5.68e-03 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
In [ ]:
P(h=0.1, rhs='-1.0-3*y*x+y*y', g='y', order=1, plot='f3.png')
P(h=0.1, rhs='-1.0-3*y*x+y*y', g='y', order=1, solver ='scimba')
P(h=0.1, rhs='-1.0-4*y*x+y*y', g='x', order=1, plot='f4.png')
P(h=0.1, rhs='-1.0-4*y*x+y*y', g='x', order=1, solver ='scimba')
P(h=0.05, rhs='1', g='0', order=1, plot='f5.png')
P(h=0.05, rhs='1 + x-x', g='0', order=1, solver ='scimba')
Solving the laplacian problem for hsize = 0.1...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo 0.1
+----------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+----------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-5 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 1.018751e-01 | | | |
| | | | h_max : 1.168628e-01 | | | |
| | | | h_min : 8.366122e-02 | | | |
| | | | n_elements : 248 | | | |
| | | | n_faces : 392 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 145 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+================+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1,0,0,1} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | -1.0+y^2-3*x*y | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-5 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 1.018751e-01 | | | | |
| | | | | h_max : 1.168628e-01 | | | | |
| | | | | h_min : 8.366122e-02 | | | | |
| | | | | n_elements : 248 | | | | |
| | | | | n_faces : 392 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 145 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : y | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-5 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 145 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Fields | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-4 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | +--------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-4 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | +--------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-4 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
+----------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 2.56e+00 epoch 0: best loss = 2.56e+00 epoch 1: best loss = 2.08e+00 epoch 2: best loss = 1.61e+00 epoch 3: best loss = 1.06e+00 epoch 4: best loss = 6.13e-01 epoch 5: best loss = 4.40e-01 epoch 6: best loss = 4.23e-01 epoch 14: best loss = 4.17e-01 epoch 15: best loss = 4.13e-01 epoch 18: best loss = 4.12e-01 epoch 19: best loss = 3.81e-01 epoch 20: best loss = 3.68e-01 epoch 21: best loss = 3.63e-01 epoch 22: best loss = 3.57e-01 epoch 23: best loss = 3.24e-01 epoch 24: best loss = 3.03e-01 epoch 25: best loss = 2.86e-01 epoch 26: best loss = 2.62e-01 epoch 27: best loss = 2.04e-01 epoch 28: best loss = 1.90e-01 epoch 31: best loss = 1.80e-01 epoch 32: best loss = 1.78e-01 epoch 33: best loss = 1.51e-01 epoch 35: best loss = 1.38e-01 epoch 37: best loss = 1.35e-01 epoch 39: best loss = 1.28e-01 epoch 40: best loss = 1.25e-01 epoch 41: best loss = 1.21e-01 epoch 42: best loss = 1.21e-01 epoch 43: best loss = 1.09e-01 epoch 46: best loss = 1.03e-01 epoch 50: best loss = 8.66e-02 epoch 51: best loss = 8.16e-02 epoch 52: best loss = 7.88e-02 epoch 57: best loss = 7.42e-02 epoch 59: best loss = 7.40e-02 epoch 60: best loss = 6.61e-02 epoch 61: best loss = 6.30e-02 epoch 87: best loss = 5.96e-02 epoch 90: best loss = 5.11e-02 epoch 104: best loss = 4.95e-02 epoch 113: best loss = 4.74e-02 epoch 126: best loss = 4.66e-02 epoch 128: best loss = 4.30e-02 epoch 134: best loss = 4.05e-02 epoch 138: best loss = 3.99e-02 epoch 141: best loss = 3.44e-02 epoch 212: best loss = 3.21e-02 epoch 225: best loss = 3.08e-02 epoch 231: best loss = 2.81e-02 epoch 233: best loss = 2.74e-02 epoch 234: best loss = 2.65e-02 epoch 238: best loss = 2.39e-02 epoch 255: best loss = 2.22e-02 epoch 258: best loss = 2.07e-02 epoch 263: best loss = 2.02e-02 epoch 273: best loss = 1.93e-02 epoch 278: best loss = 1.85e-02 epoch 281: best loss = 1.80e-02 epoch 283: best loss = 1.66e-02 epoch 297: best loss = 1.56e-02 epoch 300: best loss = 1.56e-02 epoch 306: best loss = 1.35e-02 epoch 314: best loss = 1.26e-02 epoch 327: best loss = 1.17e-02 epoch 384: best loss = 1.13e-02 epoch 395: best loss = 1.12e-02 epoch 397: best loss = 1.02e-02 epoch 403: best loss = 9.66e-03 epoch 413: best loss = 9.36e-03 epoch 417: best loss = 8.44e-03 epoch 421: best loss = 8.40e-03 epoch 426: best loss = 8.02e-03 epoch 429: best loss = 7.31e-03 epoch 448: best loss = 7.16e-03 epoch 452: best loss = 6.64e-03 epoch 462: best loss = 6.14e-03 epoch 463: best loss = 5.10e-03 epoch 484: best loss = 4.94e-03 epoch 497: best loss = 4.70e-03 epoch 500: current loss = 5.59e-03 epoch 529: best loss = 3.61e-03 epoch 542: best loss = 3.55e-03 epoch 545: best loss = 2.91e-03 epoch 552: best loss = 2.80e-03 epoch 569: best loss = 2.04e-03 epoch 599: current loss = 5.00e-03 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
Solving the laplacian problem for hsize = 0.1...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo 0.1
+----------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+----------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-7 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 1.018751e-01 | | | |
| | | | h_max : 1.168628e-01 | | | |
| | | | h_min : 8.366122e-02 | | | |
| | | | n_elements : 248 | | | |
| | | | n_faces : 392 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 145 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+================+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1,0,0,1} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | -1.0+y^2-4*y*x | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+----------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +----------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +--------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-7 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 1.018751e-01 | | | | |
| | | | | h_max : 1.168628e-01 | | | | |
| | | | | h_min : 8.366122e-02 | | | | |
| | | | | n_elements : 248 | | | | |
| | | | | n_faces : 392 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 145 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : x | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-7 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 145 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Fields | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-6 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | +--------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-6 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | +--------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-6 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| +------------------------------------------------------------------------------------------------------------------------------------------------+ |
+----------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 3.45e+00 epoch 0: best loss = 3.45e+00 epoch 1: best loss = 2.86e+00 epoch 2: best loss = 2.15e+00 epoch 3: best loss = 1.57e+00 epoch 4: best loss = 9.26e-01 epoch 5: best loss = 7.38e-01 epoch 8: best loss = 7.09e-01 epoch 9: best loss = 7.01e-01 epoch 15: best loss = 7.00e-01 epoch 16: best loss = 6.51e-01 epoch 17: best loss = 5.21e-01 epoch 19: best loss = 5.03e-01 epoch 20: best loss = 4.80e-01 epoch 21: best loss = 4.15e-01 epoch 23: best loss = 3.07e-01 epoch 24: best loss = 2.95e-01 epoch 25: best loss = 2.75e-01 epoch 27: best loss = 2.48e-01 epoch 30: best loss = 2.43e-01 epoch 31: best loss = 2.40e-01 epoch 32: best loss = 2.22e-01 epoch 35: best loss = 2.11e-01 epoch 36: best loss = 2.08e-01 epoch 39: best loss = 1.79e-01 epoch 41: best loss = 1.78e-01 epoch 43: best loss = 1.66e-01 epoch 44: best loss = 1.58e-01 epoch 47: best loss = 1.56e-01 epoch 49: best loss = 1.47e-01 epoch 52: best loss = 1.45e-01 epoch 55: best loss = 1.31e-01 epoch 57: best loss = 1.27e-01 epoch 61: best loss = 1.23e-01 epoch 64: best loss = 1.23e-01 epoch 65: best loss = 1.18e-01 epoch 69: best loss = 1.10e-01 epoch 70: best loss = 1.04e-01 epoch 74: best loss = 9.91e-02 epoch 77: best loss = 9.11e-02 epoch 78: best loss = 8.65e-02 epoch 97: best loss = 8.47e-02 epoch 100: best loss = 7.48e-02 epoch 104: best loss = 7.01e-02 epoch 105: best loss = 6.98e-02 epoch 112: best loss = 6.78e-02 epoch 115: best loss = 6.75e-02 epoch 121: best loss = 6.41e-02 epoch 123: best loss = 6.10e-02 epoch 128: best loss = 5.30e-02 epoch 129: best loss = 4.87e-02 epoch 133: best loss = 4.13e-02 epoch 141: best loss = 3.96e-02 epoch 143: best loss = 3.96e-02 epoch 144: best loss = 3.82e-02 epoch 151: best loss = 3.80e-02 epoch 153: best loss = 3.75e-02 epoch 154: best loss = 3.73e-02 epoch 158: best loss = 2.65e-02 epoch 167: best loss = 2.62e-02 epoch 169: best loss = 1.98e-02 epoch 183: best loss = 1.92e-02 epoch 197: best loss = 1.47e-02 epoch 228: best loss = 1.36e-02 epoch 234: best loss = 1.28e-02 epoch 240: best loss = 1.23e-02 epoch 259: best loss = 1.06e-02 epoch 287: best loss = 9.91e-03 epoch 297: best loss = 8.52e-03 epoch 329: best loss = 8.43e-03 epoch 330: best loss = 8.08e-03 epoch 347: best loss = 7.80e-03 epoch 362: best loss = 5.14e-03 epoch 456: best loss = 5.03e-03 epoch 479: best loss = 4.77e-03 epoch 500: current loss = 9.36e-03 epoch 502: best loss = 4.39e-03 epoch 506: best loss = 4.34e-03 epoch 510: best loss = 4.11e-03 epoch 544: best loss = 3.75e-03 epoch 599: current loss = 7.38e-03 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
Solving the laplacian problem for hsize = 0.05...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo 0.05
+------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-9 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 5.095217e-02 | | | |
| | | | h_max : 6.887751e-02 | | | |
| | | | h_min : 4.286120e-02 | | | |
| | | | n_elements : 952 | | | |
| | | | n_faces : 1468 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 517 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +--------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +--------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +----------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +----------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +----------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +----------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+============+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1,0,0,1} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | 1 | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +----------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +--------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +----------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-9 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 5.095217e-02 | | | | |
| | | | | h_max : 6.887751e-02 | | | | |
| | | | | h_min : 4.286120e-02 | | | | |
| | | | | n_elements : 952 | | | | |
| | | | | n_faces : 1468 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 517 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : 0 | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-9 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 517 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Fields | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-8 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | +--------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-8 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | +--------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +--------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-8 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| +--------------------------------------------------------------------------------------------------------------------------------------------+ |
+------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 8.02e-01 epoch 0: best loss = 8.02e-01 epoch 1: best loss = 5.90e-01 epoch 2: best loss = 3.73e-01 epoch 3: best loss = 1.55e-01 epoch 4: best loss = 1.12e-01 epoch 8: best loss = 9.63e-02 epoch 13: best loss = 9.29e-02 epoch 19: best loss = 9.24e-02 epoch 28: best loss = 9.14e-02 epoch 33: best loss = 9.03e-02 epoch 40: best loss = 8.95e-02 epoch 41: best loss = 8.90e-02 epoch 56: best loss = 8.86e-02 epoch 63: best loss = 8.79e-02 epoch 99: best loss = 8.74e-02 epoch 109: best loss = 8.61e-02 epoch 113: best loss = 8.26e-02 epoch 116: best loss = 8.22e-02 epoch 118: best loss = 7.91e-02 epoch 119: best loss = 7.69e-02 epoch 120: best loss = 6.89e-02 epoch 121: best loss = 6.53e-02 epoch 122: best loss = 5.95e-02 epoch 123: best loss = 5.31e-02 epoch 124: best loss = 4.75e-02 epoch 127: best loss = 4.74e-02 epoch 128: best loss = 3.89e-02 epoch 129: best loss = 3.27e-02 epoch 130: best loss = 2.90e-02 epoch 140: best loss = 2.89e-02 epoch 141: best loss = 2.45e-02 epoch 142: best loss = 2.41e-02 epoch 147: best loss = 2.18e-02 epoch 151: best loss = 2.16e-02 epoch 155: best loss = 2.14e-02 epoch 156: best loss = 1.94e-02 epoch 157: best loss = 1.87e-02 epoch 161: best loss = 1.81e-02 epoch 163: best loss = 1.78e-02 epoch 164: best loss = 1.76e-02 epoch 165: best loss = 1.70e-02 epoch 166: best loss = 1.65e-02 epoch 167: best loss = 1.62e-02 epoch 170: best loss = 1.53e-02 epoch 171: best loss = 1.40e-02 epoch 174: best loss = 1.33e-02 epoch 175: best loss = 1.25e-02 epoch 177: best loss = 1.19e-02 epoch 178: best loss = 1.11e-02 epoch 181: best loss = 1.03e-02 epoch 182: best loss = 9.72e-03 epoch 183: best loss = 9.64e-03 epoch 185: best loss = 8.48e-03 epoch 186: best loss = 8.00e-03 epoch 192: best loss = 7.37e-03 epoch 194: best loss = 6.94e-03 epoch 196: best loss = 6.29e-03 epoch 197: best loss = 6.17e-03 epoch 202: best loss = 6.13e-03 epoch 207: best loss = 4.96e-03 epoch 217: best loss = 4.39e-03 epoch 231: best loss = 4.38e-03 epoch 236: best loss = 4.06e-03 epoch 242: best loss = 4.00e-03 epoch 243: best loss = 3.97e-03 epoch 245: best loss = 3.60e-03 epoch 258: best loss = 3.30e-03 epoch 259: best loss = 3.08e-03 epoch 266: best loss = 2.78e-03 epoch 288: best loss = 2.59e-03 epoch 294: best loss = 2.58e-03 epoch 325: best loss = 2.53e-03 epoch 329: best loss = 2.43e-03 epoch 337: best loss = 2.31e-03 epoch 338: best loss = 2.21e-03 epoch 339: best loss = 2.06e-03 epoch 340: best loss = 2.06e-03 epoch 342: best loss = 1.89e-03 epoch 352: best loss = 1.75e-03 epoch 370: best loss = 1.72e-03 epoch 376: best loss = 1.64e-03 epoch 383: best loss = 1.63e-03 epoch 385: best loss = 1.41e-03 epoch 398: best loss = 1.38e-03 epoch 405: best loss = 1.19e-03 epoch 500: current loss = 1.52e-03 epoch 511: best loss = 8.45e-04 epoch 534: best loss = 7.84e-04 epoch 599: current loss = 1.08e-03 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
In [ ]:
# # 2D with varying anisotropy
#P = Poisson(dim = 2)
P(h=0.1, diff='{1.0,0,0,x*y}', rhs='1', plot='d1.png')
P(h=0.1, diff='(1+x,0,0,1+y)', rhs='1+ x-x', solver='scimba')
P(h=0.1, diff='{x,y,-y,x+y}', rhs='1', plot='d3.png')
P(h=0.1, diff='(x,y,-y,x+y)', rhs='1+ x-x', solver='scimba')
Solving the laplacian problem for hsize = 0.1...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo 0.1
+---------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+---------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-11 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 1.018751e-01 | | | |
| | | | h_max : 1.168628e-01 | | | |
| | | | h_min : 8.366122e-02 | | | |
| | | | n_elements : 248 | | | |
| | | | n_faces : 392 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 145 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +-----------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +-----------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +-------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +-------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +---------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +---------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +-----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +-----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +-------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +-------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+---------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+===============+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1.0,0,0,x*y} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+---------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | 1 | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+---------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +-------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +-----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +---------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +-------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-11 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 1.018751e-01 | | | | |
| | | | | h_max : 1.168628e-01 | | | | |
| | | | | h_min : 8.366122e-02 | | | | |
| | | | | n_elements : 248 | | | | |
| | | | | n_faces : 392 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 145 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : 0 | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-11 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 145 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +-------------------------------------------------+ | |
| | | Fields | | |
| | +-------------------------------------------------+ | |
| | | +---------------------------------------------+ | | |
| | | | Potential | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-10 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-10 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-10 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | +-------------------------------------------------+ | |
| +-----------------------------------------------------------------------------------------------------------------------------------------------+ |
+---------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 1.00e+00 epoch 0: best loss = 1.00e+00 epoch 1: best loss = 6.21e-01 epoch 2: best loss = 2.75e-01 epoch 3: best loss = 9.68e-02 epoch 13: best loss = 9.06e-02 epoch 19: best loss = 8.91e-02 epoch 35: best loss = 8.70e-02 epoch 40: best loss = 8.59e-02 epoch 56: best loss = 8.59e-02 epoch 86: best loss = 8.55e-02 epoch 90: best loss = 8.35e-02 epoch 93: best loss = 8.00e-02 epoch 96: best loss = 7.82e-02 epoch 97: best loss = 7.46e-02 epoch 98: best loss = 7.09e-02 epoch 99: best loss = 6.77e-02 epoch 100: best loss = 6.40e-02 epoch 101: best loss = 6.28e-02 epoch 102: best loss = 5.69e-02 epoch 116: best loss = 5.50e-02 epoch 117: best loss = 5.49e-02 epoch 118: best loss = 5.28e-02 epoch 120: best loss = 5.27e-02 epoch 121: best loss = 5.24e-02 epoch 122: best loss = 5.11e-02 epoch 123: best loss = 5.08e-02 epoch 124: best loss = 4.94e-02 epoch 125: best loss = 4.87e-02 epoch 128: best loss = 4.66e-02 epoch 131: best loss = 4.66e-02 epoch 132: best loss = 4.49e-02 epoch 134: best loss = 4.14e-02 epoch 136: best loss = 3.99e-02 epoch 137: best loss = 3.83e-02 epoch 138: best loss = 3.47e-02 epoch 139: best loss = 3.16e-02 epoch 140: best loss = 2.93e-02 epoch 141: best loss = 2.58e-02 epoch 142: best loss = 2.48e-02 epoch 145: best loss = 2.48e-02 epoch 148: best loss = 2.42e-02 epoch 149: best loss = 2.34e-02 epoch 150: best loss = 2.21e-02 epoch 151: best loss = 2.01e-02 epoch 158: best loss = 1.97e-02 epoch 159: best loss = 1.83e-02 epoch 161: best loss = 1.63e-02 epoch 164: best loss = 1.62e-02 epoch 165: best loss = 1.61e-02 epoch 167: best loss = 1.60e-02 epoch 168: best loss = 1.34e-02 epoch 175: best loss = 1.26e-02 epoch 176: best loss = 1.23e-02 epoch 177: best loss = 1.20e-02 epoch 179: best loss = 1.18e-02 epoch 180: best loss = 1.12e-02 epoch 181: best loss = 1.04e-02 epoch 182: best loss = 9.99e-03 epoch 184: best loss = 8.24e-03 epoch 185: best loss = 8.17e-03 epoch 186: best loss = 8.16e-03 epoch 187: best loss = 6.54e-03 epoch 191: best loss = 6.41e-03 epoch 192: best loss = 5.70e-03 epoch 194: best loss = 5.60e-03 epoch 195: best loss = 5.42e-03 epoch 197: best loss = 4.30e-03 epoch 198: best loss = 3.92e-03 epoch 201: best loss = 3.88e-03 epoch 202: best loss = 3.58e-03 epoch 204: best loss = 3.23e-03 epoch 210: best loss = 2.83e-03 epoch 216: best loss = 2.56e-03 epoch 218: best loss = 2.52e-03 epoch 221: best loss = 2.32e-03 epoch 228: best loss = 2.30e-03 epoch 230: best loss = 2.24e-03 epoch 234: best loss = 2.21e-03 epoch 236: best loss = 1.92e-03 epoch 240: best loss = 1.79e-03 epoch 242: best loss = 1.35e-03 epoch 265: best loss = 1.24e-03 epoch 267: best loss = 1.13e-03 epoch 278: best loss = 8.51e-04 epoch 347: best loss = 6.39e-04 epoch 386: best loss = 5.54e-04 epoch 415: best loss = 5.41e-04 epoch 469: best loss = 3.79e-04 epoch 500: current loss = 8.55e-04 epoch 508: best loss = 3.31e-04 epoch 520: best loss = 3.30e-04 epoch 599: current loss = 4.08e-04 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
Solving the laplacian problem for hsize = 0.1...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/omega-2.geo 0.1
+--------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-2d-p1 - Use Case Study |
+--------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-2d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-2d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-13 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 2 | | | |
| | | | h_average : 1.018751e-01 | | | |
| | | | h_max : 1.168628e-01 | | | |
| | | | h_min : 8.366122e-02 | | | |
| | | | n_elements : 248 | | | |
| | | | n_faces : 392 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 145 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 2 | | | |
| | | | shape : Simplex_2_1_2 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +----------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +----------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +--------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +--------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+--------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+==============+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {x,y,-y,x+y} | physics_poisson_eq_poisson_eq_c | tensor2 [2x2] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+--------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | 1 | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+--------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +----------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +--------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-13 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 2 | | | | |
| | | | | h_average : 1.018751e-01 | | | | |
| | | | | h_max : 1.168628e-01 | | | | |
| | | | | h_min : 8.366122e-02 | | | | |
| | | | | n_elements : 248 | | | | |
| | | | | n_faces : 392 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 145 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 2 | | | | |
| | | | | shape : Simplex_2_1_2 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : 0 | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-13 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 3 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +--------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +--------------+ | | | |
| | | | | nDof : 145 | | | | |
| | | | +--------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +-------------------------------------------------+ | |
| | | Fields | | |
| | +-------------------------------------------------+ | |
| | | +---------------------------------------------+ | | |
| | | | Potential | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-12 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-12 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-12 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | +-------------------------------------------------+ | |
| +----------------------------------------------------------------------------------------------------------------------------------------------+ |
+--------------------------------------------------------------------------------------------------------------------------------------------------+
Solving using Scimba >> load network /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth network was not loaded from file: training needed epoch 0: current loss = 9.97e-01 epoch 0: best loss = 9.97e-01 epoch 1: best loss = 6.29e-01 epoch 2: best loss = 3.18e-01 epoch 3: best loss = 1.26e-01 epoch 7: best loss = 1.20e-01 epoch 8: best loss = 9.27e-02 epoch 13: best loss = 8.91e-02 epoch 27: best loss = 8.87e-02 epoch 41: best loss = 8.56e-02 epoch 55: best loss = 8.49e-02 epoch 56: best loss = 8.41e-02 epoch 58: best loss = 8.23e-02 epoch 61: best loss = 7.83e-02 epoch 65: best loss = 7.65e-02 epoch 66: best loss = 7.53e-02 epoch 68: best loss = 7.03e-02 epoch 69: best loss = 6.80e-02 epoch 70: best loss = 6.10e-02 epoch 71: best loss = 5.86e-02 epoch 72: best loss = 5.15e-02 epoch 73: best loss = 4.65e-02 epoch 74: best loss = 4.56e-02 epoch 78: best loss = 4.17e-02 epoch 79: best loss = 3.38e-02 epoch 80: best loss = 3.37e-02 epoch 81: best loss = 3.08e-02 epoch 84: best loss = 3.04e-02 epoch 94: best loss = 2.87e-02 epoch 97: best loss = 2.81e-02 epoch 99: best loss = 2.69e-02 epoch 109: best loss = 2.67e-02 epoch 112: best loss = 2.53e-02 epoch 115: best loss = 2.49e-02 epoch 117: best loss = 2.25e-02 epoch 120: best loss = 2.18e-02 epoch 124: best loss = 2.00e-02 epoch 132: best loss = 1.96e-02 epoch 134: best loss = 1.94e-02 epoch 139: best loss = 1.93e-02 epoch 154: best loss = 1.86e-02 epoch 155: best loss = 1.82e-02 epoch 158: best loss = 1.80e-02 epoch 162: best loss = 1.75e-02 epoch 164: best loss = 1.74e-02 epoch 167: best loss = 1.73e-02 epoch 168: best loss = 1.72e-02 epoch 170: best loss = 1.68e-02 epoch 173: best loss = 1.66e-02 epoch 175: best loss = 1.65e-02 epoch 176: best loss = 1.61e-02 epoch 183: best loss = 1.49e-02 epoch 184: best loss = 1.40e-02 epoch 188: best loss = 1.39e-02 epoch 193: best loss = 1.25e-02 epoch 197: best loss = 1.18e-02 epoch 201: best loss = 1.03e-02 epoch 203: best loss = 1.03e-02 epoch 205: best loss = 1.02e-02 epoch 206: best loss = 9.01e-03 epoch 207: best loss = 7.78e-03 epoch 208: best loss = 7.27e-03 epoch 209: best loss = 6.71e-03 epoch 210: best loss = 6.14e-03 epoch 211: best loss = 5.98e-03 epoch 213: best loss = 5.18e-03 epoch 223: best loss = 5.08e-03 epoch 224: best loss = 4.57e-03 epoch 228: best loss = 4.51e-03 epoch 233: best loss = 3.94e-03 epoch 234: best loss = 3.82e-03 epoch 238: best loss = 3.27e-03 epoch 239: best loss = 3.17e-03 epoch 242: best loss = 2.68e-03 epoch 262: best loss = 2.62e-03 epoch 296: best loss = 2.44e-03 epoch 298: best loss = 2.44e-03 epoch 303: best loss = 2.25e-03 epoch 304: best loss = 2.08e-03 epoch 322: best loss = 1.92e-03 epoch 323: best loss = 1.76e-03 epoch 337: best loss = 1.40e-03 epoch 377: best loss = 1.32e-03 epoch 424: best loss = 1.12e-03 epoch 452: best loss = 1.07e-03 epoch 480: best loss = 9.04e-04 epoch 494: best loss = 8.88e-04 epoch 500: current loss = 1.25e-03 epoch 502: best loss = 8.40e-04 epoch 524: best loss = 5.53e-04 epoch 562: best loss = 5.45e-04 epoch 584: best loss = 4.28e-04 epoch 599: current loss = 7.57e-04 load network: /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1/networks/test.pth jsbdbshd network loaded
In [ ]:
# 3D with varying anisotropy and non-homogeneous Dirichlet BC
P = Poisson(dim = 3)
P(h=0.08, diff='{1,0,0,0,x+1,0,0,0,1+x*y}', g = 'x', rhs='x*y*z', geofile = 'geo/cube.geo', plot='3d.png')
Solving the laplacian problem for hsize = 0.08...
[loadMesh] Loading mesh in format geo+msh: "/workspaces/2024-m1-scimba-feelpp/geo/cube.geo"
[loadMesh] Use default geo desc: /workspaces/2024-m1-scimba-feelpp/geo/cube.geo 0.08
+---------------------------------------------------------------------------------------------------------------------------------------------------------------+
| Toolbox::cfpdes-3d-p1 - Use Case Study |
+---------------------------------------------------------------------------------------------------------------------------------------------------------------+
| +--------------------------------------------------------------------------------------+ |
| | Environment | |
| +--------------------------------------------------------------------------------------+ |
| | prefix : cfpdes | |
| | keyword : cfpdes-3d-p1 | |
| | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | |
| | number of processus : 1 | |
| +--------------------------------------------------------------------------------------+ |
| +---------------------------+ |
| | Materials Properties | |
| +---------------------------+ |
| | number of materials : 1 | |
| +---------------------------+ |
| +------------------------------------------------+ |
| | Meshes | |
| +------------------------------------------------+ |
| | +--------------------------------------------+ | |
| | | Mesh : cfpdes-3d-p1 | | |
| | +--------------------------------------------+ | |
| | | Discretization : /Mesh/object-15 | | |
| | | +----------------------------------------+ | | |
| | | | Import configuration | | | |
| | | +----------------------------------------+ | | |
| | | | geo-filename : | | | |
| | | | hsize : 1.000000e-01 | | | |
| | | | generate-partitioning : 0 | | | |
| | | +----------------------------------------+ | | |
| | | +-------------------------------+ | | |
| | | | Discretization | | | |
| | | +-------------------------------+ | | |
| | | | dim : 3 | | | |
| | | | h_average : 1.254175e-01 | | | |
| | | | h_max : 1.686836e-01 | | | |
| | | | h_min : 7.692308e-02 | | | |
| | | | n_edges : 13802 | | | |
| | | | n_elements : 10288 | | | |
| | | | n_faces : 21787 | | | |
| | | | n_partition : 1 | | | |
| | | | n_points : 2304 | | | |
| | | | order : 1 | | | |
| | | | real_dim : 3 | | | |
| | | | shape : Simplex_3_1_3 | | | |
| | | +-------------------------------+ | | |
| | +--------------------------------------------+ | |
| +------------------------------------------------+ |
| +----------------------------------+ |
| | Algebraic Solver | |
| +----------------------------------+ |
| | +-------------------+ | |
| | | Backend | | |
| | +-------------------+ | |
| | | prefix : cfpdes | | |
| | | type : petsc | | |
| | +-------------------+ | |
| | +-----------------------------+ | |
| | | KSP | | |
| | +-----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | dtol : 1.000000e+05 | | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | type : gmres | | |
| | +-----------------------------+ | |
| | +----------------------------+ | |
| | | SNES | | |
| | +----------------------------+ | |
| | | atol : 1.000000e-50 | | |
| | | maxit : 50 | | |
| | | reuse-jac : 0 | | |
| | | rtol : 1.000000e-08 | | |
| | | stol : 1.000000e-08 | | |
| | +----------------------------+ | |
| | +-----------------------------+ | |
| | | KSP in SNES | | |
| | +-----------------------------+ | |
| | | maxit : 1000 | | |
| | | reuse-prec : 0 | | |
| | | rtol : 1.000000e-05 | | |
| | +-----------------------------+ | |
| | +------------------------------+ | |
| | | PC | | |
| | +------------------------------+ | |
| | | mat-solver-package : mumps | | |
| | | type : lu | | |
| | +------------------------------+ | |
| +----------------------------------+ |
| +-----------------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | Toolbox Coefficient Form PDE : poisson_eq | |
| +-----------------------------------------------------------------------------------------------------------------------------------------------------------+ |
| | +--------------------------------------------------------------------------------------+ | |
| | | Environment | | |
| | +--------------------------------------------------------------------------------------+ | |
| | | prefix : cfpdes.poisson_eq | | |
| | | keyword : poisson_eq | | |
| | | root repository : /workspaces/2024-m1-scimba-feelpp/feelppdb/feelpp_cfpde/np_1 | | |
| | | number of processus : 1 | | |
| | +--------------------------------------------------------------------------------------+ | |
| | +-------------------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | Physics | | |
| | +-------------------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | | +---------------------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | poisson_eq | | | |
| | | +---------------------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | | | +-----------------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | poisson_eq | | | | |
| | | | +-----------------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | | | modeling : GenericPDE | | | | |
| | | | | type : poisson_eq | | | | |
| | | | | name : poisson_eq | | | | |
| | | | | +-------------------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | Parameters | | | | | |
| | | | | +-------------------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | | | +------+---------------------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | Name | Expression | Symbol | Shape | Components | | | | | | |
| | | | | | +======+===========================+=================================+===============+==================================================+ | | | | | |
| | | | | | | c | {1,0,0,0,1+x,0,0,0,1+x*y} | physics_poisson_eq_poisson_eq_c | tensor2 [3x3] | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | Symbol | Indices | | | | | | | |
| | | | | | | | | | | +====================================+=========+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_00 | 0,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_01 | 0,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_02 | 0,2 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_10 | 1,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_11 | 1,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_12 | 1,2 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_20 | 2,0 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_21 | 2,1 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | | | | | | | physics_poisson_eq_poisson_eq_c_22 | 2,2 | | | | | | | |
| | | | | | | | | | | +------------------------------------+---------+ | | | | | | |
| | | | | | +------+---------------------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | | | f | y*z*x | physics_poisson_eq_poisson_eq_f | scalar | | | | | | | |
| | | | | | +------+---------------------------+---------------------------------+---------------+--------------------------------------------------+ | | | | | |
| | | | | +-------------------------------------------------------------------------------------------------------------------------------------------+ | | | | |
| | | | +-----------------------------------------------------------------------------------------------------------------------------------------------+ | | | |
| | | +---------------------------------------------------------------------------------------------------------------------------------------------------+ | | |
| | +-------------------------------------------------------------------------------------------------------------------------------------------------------+ | |
| | +------------------------------------------------+ | |
| | | Meshes | | |
| | +------------------------------------------------+ | |
| | | +--------------------------------------------+ | | |
| | | | Mesh : poisson_eq | | | |
| | | +--------------------------------------------+ | | |
| | | | Discretization : /Mesh/object-15 | | | |
| | | | +----------------------------------------+ | | | |
| | | | | Import configuration | | | | |
| | | | +----------------------------------------+ | | | |
| | | | | geo-filename : | | | | |
| | | | | hsize : 1.000000e-01 | | | | |
| | | | | generate-partitioning : 0 | | | | |
| | | | +----------------------------------------+ | | | |
| | | | +-------------------------------+ | | | |
| | | | | Discretization | | | | |
| | | | +-------------------------------+ | | | |
| | | | | dim : 3 | | | | |
| | | | | h_average : 1.254175e-01 | | | | |
| | | | | h_max : 1.686836e-01 | | | | |
| | | | | h_min : 7.692308e-02 | | | | |
| | | | | n_edges : 13802 | | | | |
| | | | | n_elements : 10288 | | | | |
| | | | | n_faces : 21787 | | | | |
| | | | | n_partition : 1 | | | | |
| | | | | n_points : 2304 | | | | |
| | | | | order : 1 | | | | |
| | | | | real_dim : 3 | | | | |
| | | | | shape : Simplex_3_1_3 | | | | |
| | | | +-------------------------------+ | | | |
| | | +--------------------------------------------+ | | |
| | +------------------------------------------------+ | |
| | +---------------------------------+ | |
| | | Boundary Conditions | | |
| | +---------------------------------+ | |
| | | +-----------------------------+ | | |
| | | | Dirichlet | | | |
| | | +-----------------------------+ | | |
| | | | +-------------------------+ | | | |
| | | | | g | | | | |
| | | | +-------------------------+ | | | |
| | | | | method : elimination | | | | |
| | | | | expr : x | | | | |
| | | | | markers : Gamma_D | | | | |
| | | | +-------------------------+ | | | |
| | | +-----------------------------+ | | |
| | +---------------------------------+ | |
| | +------------------------------------+ | |
| | | Function Spaces | | |
| | +------------------------------------+ | |
| | | +--------------------------------+ | | |
| | | | Potential | | | |
| | | +--------------------------------+ | | |
| | | | mesh : /Mesh/object-15 | | | |
| | | | nSpace : 1 | | | |
| | | | +----------------------------+ | | | |
| | | | | Basis | | | | |
| | | | +----------------------------+ | | | |
| | | | | is_continuous : 1 | | | | |
| | | | | nComponents : 1 | | | | |
| | | | | nComponents1 : 1 | | | | |
| | | | | nComponents2 : 1 | | | | |
| | | | | nLocalDof : 4 | | | | |
| | | | | name : lagrange | | | | |
| | | | | order : 1 | | | | |
| | | | | shape : scalar | | | | |
| | | | +----------------------------+ | | | |
| | | | +---------------+ | | | |
| | | | | Dof Table | | | | |
| | | | +---------------+ | | | |
| | | | | nDof : 2304 | | | | |
| | | | +---------------+ | | | |
| | | +--------------------------------+ | | |
| | +------------------------------------+ | |
| | +-------------------------------------------------+ | |
| | | Fields | | |
| | +-------------------------------------------------+ | |
| | | +---------------------------------------------+ | | |
| | | | Potential | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u | | | |
| | | | function space : /FunctionSpace/object-14 | | | |
| | | | name : Potential | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_previous | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_previous | | | |
| | | | function space : /FunctionSpace/object-14 | | | |
| | | | name : Potential_previous | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | | +---------------------------------------------+ | | |
| | | | Potential_remove_trial | | | |
| | | +---------------------------------------------+ | | |
| | | | base symbol : u_rt | | | |
| | | | function space : /FunctionSpace/object-14 | | | |
| | | | name : Potential_remove_trial | | | |
| | | | prefix symbol : poisson_eq | | | |
| | | +---------------------------------------------+ | | |
| | +-------------------------------------------------+ | |
| +-----------------------------------------------------------------------------------------------------------------------------------------------------------+ |
+---------------------------------------------------------------------------------------------------------------------------------------------------------------+